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Abstract. The spectral problem for matrices with a block- hierarchical structure is often 
considered in context of the theory of complex systems. In the present article, a new class 
of matrices with a block-rectangular non-symmetric hierarchical structure is introduced and 
the corresponding spectral problem is investigated. Using these results we study a model of 
error generation in information sequence where such block-rectangular hierarchical matrices 
appear in a natural way. 



1. Introduction 

Block-hierarchical matrices are distingushed by their specific structure of hierarchically 
nested growing blocks placed along the diagonal, where each (sub-)block is again a block- 
hierarchical matrix itself, see Ex.1. Essential necessity in mathematical objects of such 
a nature arose from studying of spin-glass models at the end of the 70's [1,2]. In the 
most studied Sherrington-Kirkpatrick model [3] the block-hierarchical matrices appear as 
solutions of a saddle-point equation in the calculation of the free energy function by means 
of the replica method. The entries of these matrices are interpreted as the overlap parameter 
between different states of the spin-glass system which correspond to the minima of the free 
energy function [2]. Since the set of these minima has a hierarchically nested structure the 
notion of ultrametric space [4] arises in a natural way. This, in turn, is reflected in the 
block-hierarchical structure of the correlation matrices. 
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Example 1 A block-hierarchical matrix for p = 2: 
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The 2x2 blocks placed along the diagonal are embedded into 4x4 blocks and so on. The blocks of 
the level-number 7, 7 = 1, 2, . . . have the size 2 7 x 2 7 . 

The distance, d(x,y), denned on the "ordinary" metric spaces must satisfy the triangle 
inequality: d(x, y) < d(x, z) + d(z, y) where x, y and z are elements of the space. For spaces 
with ultrametric distance d u (x,y) some stronger relation - the strong triangle inequality - 
takes place: d u (x,y) < max {d u (x, z), d u (z, y)}. The metric structure of these spaces can be 
conveniently described by means of a directed distances graph having a tree-like structure. 
Elements of the ultrametric space are associated with the endpoints of the graph. The 
distance between two elements x, y is determined by the hierarchical level j(x,y) of their 
first common ancestor on the tree, see Fig.l. The ultrametric distance d u (x,y) is a positive 
growing function of 7 with a supplementary constraint: d u (x, y) — if and only if x — y. 

Beyond the theory of spin glasses a number of works have been devoted to the 
investigation of non-equilibrium dynamics of various complex systems [5,6] with ultrametric 
structure of the free energy landscape. In these studies the block-hierarchical matrices 
typically describe the probabilities of transitions between different states of the system 
identified with the points of the ultrametric space. The first mathematical model of such 
a kind was introduced by Ogielski and Stein in [7]. They considered a stochastic process 
which can be interpreted as "ultrametric diffusion". In this model the diffusion occurs on 
the end-points of a directed regular tree with the branching parameter p (the number of 
outcoming branches at each vertex) and is governed by the following equation 
df{t) 



dt 



= Qf(t). 



(i) 



Here the i-th element of the vector f(t) determines the probability to find the system at the 
site i at a given instant of time t. The hopping matrix Q of the block-hierarchical type (see 
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Figure 1. Distance graph of an ultrametric space consisting of p r elements (depicted 
by the numerated endpoints of the graph) with branching parameter p = 2 and 
parameter r = 3. The distance between two distinguishable elements is determined 
by the value 7 - the level of hierarchy of their first common ancestor. For instance, 
distance between elements denoted by 2 and 4 is given by 7 = 2. 



Ex.1) is composed of the nested blocks of exponentially growing sizes p 1 x p 1 , 7 = 1, 2, . . . , r. 
The solution of the lattice equation (1) reduces to the spectral problem for the matrix Q. 
Ogielski and Stein succeeded to calculate eigenvalues and eigenvectors of matrices Q (for 
p = 2) in the limit r — >• 00 and found the solution of the initial-state-decay problem. Later, 
it was shown in [8, 9] that whenever the distance graph possesses a regular structure, a 
natural framework for the analysis of eq.(l) is provided by the spectral theory of p-adic 
pseudodifferential operators [10]. However, for block-hierarchical matrices of a general form, 
where the distance graph is irregularly branched, the corresponding ultrametric space does 
not fit into the structure of the field of p-adic numbers and p-adic formalism cannot directly 
be applied [11]. 

In this work our attention is focused on a new class of non-symmetric block-hierarchical 
matrices Q. Unlike the conventional block-hierarchical matrices, Q. have rectangular (rather 
then square) structure of subblocks. More precisely, p r x p r matrix Q is composed of p 
identical p r x p r " x rectangular sub-matrices O: 



The sub-matrix O has the structure similar to conventional block-hierarchical matrices with 
the elements replaced by columns of the height p, as visualised in Ex.2. Accordingly, the 
nested rectangular sub-blocks in the matrix O are of the dimensions p 7 x p 7_1 , 7 = 2, . . . , r. 




block-rectangular hierarchical matrices 



Example 2 An 8 
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block-rectangular hierarchical matrix Q with the parameters p = 2, r 
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Here q2 is the block of the size 4x2. 

In the present paper we aim to describe the spectral properties of the block-rectangular 
hierarchical matrices. As a n application, we also provide a simple error generating model, 
where such matrices appear in a natural way. In particular, the knowledge of the spectrum of 
the block-rectangular hierarchical matrices makes it possible to evaluate the average number 
of generated errors and the corresponding variance. 

The paper is organised as follows. In the next section we introduce tensor product 
representation for standard block-hierarchical matrices and recall their spectral properties. In 
the section 3 we use this representation to introduce the class of block-rectangular hierarchical 
matrices and study the corresponding spectral problem. In the section 4 the error generating 
model is considered. The concluding remarks are given in the section 5. 



2. Block-hierarchical matrices. Tensor product representation 

The original construction of block-hierarchical matrices given in the paper [1] by Parisi might 
be described in the following way: 

Take g 7 , 7 = 0, . . . , r be a sequence ofr > 1 arbitrary numbers and take n 7 be a sequence 
of positive integer numbers, such that the ratios n 7 /n 7 _i are integers for 7 > 1. Then a non- 
diagonal element Q a b, a ^b of the n r xn r block-hierarchical matrix Q is equal to g 7 whenever 
t 



a 




' b ' 
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" b ' 
























n 7 




n 7 



The diagonal elements Q aa = q are allowed to be arbitrary, 
i The brackets [•] denote the integer part, such that for all real x 
\x] - 1 < x < \x] . 
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In what follows we focus on a class of regular block-hierarchical matrices, where n 7 /n 7 _x 
is the same for all 7. Below we introduce a tensor product representation for such matrices 
and reproduce the known results on their spectrum. 

2.1. Tensor product representation 

In many applications the sequence n 7 in the above definition is chosen to be a geometric 
progression, 

n 7 =p 7 , 7 = 0, l,...,r. (2) 

Among a wide variety of block-hierarchical matrices this condition singles out those which 
have regularly growing blocks. These matrices correspond to the ultrametric spaces with 
regularly branching distances graphs. For instance, the Ogilski-Stein model of ultrametric 
diffusion [7] was formulated for p — 2, see Ex.1. In p-adic models [8] the parameter p is 
assumed to be a prime number as it is implied in the p-adic formalism [10]. In the present 
paper this condition is relaxed and p is an arbitrary integer larger then 1, while r is a fixed 
integer. 

The structure of block-hierarchical matrices makes it possible to present them in the 
form of linear combination of block-diagonal matrices. Indeed, let us introduce p r x p r block- 
diagonal matrices 5 7 , 7 = 0, 1, . . . , r. Each 5 7 is composed of p r ~ 1 identical p 7 x p 7 blocks 
with all entries equal to one (So is nothing but the identity matrix). It is easy to check, that 
a block-hierarchical matrix Q can then be constructed in the form of linear combination of 
S^ si 



r 



Q = ^a 7 5 7 . (3) 

7=0 

Each matrix S 7 , in turn, can be represented as the tensor product: 



5 7 = 1 . . I a . . . s , (4) 

r— 7 7 

composed of r — 7 identity matrices 1 and 7 projection matrices s of the same size p x p: 
(l ... l\ 



1 

s — - 

p 



1 ... 1 



(5) 



V 1 ••• 1/ 

Here the normalisation factor 1/p is added to satisfy the projection property 
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The representation (4) together with the property (6) implies that the product of 



matrices S 7 is 

SrySy = Sry'Sry = <S maX { J j j' } . (7) 

In particular, it means that both products and sums of block-hierarchical matrices are again 
block-hierarchical matrices. 

Remark It worth mentioning that due to the relation (7), matrix Q can be cast in a product and 
exponential forms: 

r r 

Q = Y^ <h s i = e E ^ ^ = JJ (1 + fr 7 S 7 ) . 

7=0 7=0 

Equating of coefficients at different in each of the relations above results in the following 
connection formulae between a 7 , 6 7 , c 7 and g 7 : 

r 

q r = a r p~ r , q 1 = ^2 a ixP~^ ■, a i = P 7 («7 ~ Qj+i), b~ 1 = e c "> -\, 7 = 0, . . . , r - 1. (8) 

7-1 /7-I 

a = 1 + b , a 7 = 6 7 (1 + 6 M ) , 6 7 = a 7 / ^ a M , 7 = 1, . . . , r - 1. 

m=o / m=o 



2.2. Spectral properties of Q 

Let us turn now our attention to the solution of the spectral problem for block-hierarchical 
matrices. To this end it is convenient to consider an operator formulation of the problem. 
To pass from matrices to operators we, first, introduce a p-dimensional Hilbert space h and 
the associated p r -dimensional space ri with the tensor product structure 

ri = h®h® ...®h. (9) 

r 

In what follows we assume that a set of the orthonormal vectors 

\ej)eh, j e {0, l,...,p- 1}; (ei\ej) = 5ij (10) 
forms a basis in h. Correspondingly, the set of the vectors 

\ e j > : = l e ii) ® \ e j 2 ) ® • • • ® |e Jr ), j'fc e {0, 1}, (11) 

indexed by j = 1 + $^7=1 hV r ~ 1 running from 1 to p r , provides an orthonormal basis on ri. 
Let us now replace the matrix s in eq.(4) with the operator s acting as a projector onto the 
one-dimensional subspace of h 

1 

*l^> = je{0,l,...,p-l}. (12) 

^ i=0 
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By this definition the matrix s, eq. (5), can be seen as the representation of s in the basis (10), 
i.e. Sji = (j\s\i). 

Having introduced the basic notations we define a family of operators V reproducing 
the block-hierarchical matrices. 

Definition 1 Let h be a p- dimensional Hilbert and let H be the associated p r -dimensional 
Hilbert space with the structure of the tensor product (9). Let s be the projection operator 
onto one- dimensional subspace of h, as defined by (12). The family V contains all linear 
operators acting on vectors u> = u)\ (8) u 2 <S> • • • <S> oj r E Ti. according to the rule 

r 

Qu> = cijSj u>, S 7 u> — uji <S> ■ ■ ■ <S> w r _ 7 ® sw r _ 7+1 (g> . . . <g> su r , (13) 

7=0 

where a 7 's are arbitrary numbers. 

It is easy to see that any Q E V takes a block-hierarchical form in the basis (11). 

The tensor product form of operators Q eV makes the solution of the spectral problem 
straightforward. Denote j — 1, . . . ,p — 1 the orthogonal eigenvectors of s, such that: 

J|0) = 1.|0>, 

s\j)=0-\j), j E {1,. ..,£> — 1}. 
Explicitly these eigenvectors can be expressed in the basis (10) as 

\j) = 4= E c ^|e fc ), cf = e 2 ^ k, j E {0, 1, . . . ,p - 1}. (14) 

(Note, that such a choice of the coefficients cf is prompted by the structure of multiplicative 
characters of the field of p-adic numbers [10].) Then from eq.(13) it follows that any tensor 
product state 

\j1j2 ■ ■ ■ 3r >:= \ji) <8> |j2> ® • • • ® \j r ), jk G {0, 1 , . . . , p - 1 } , (15) 

which is composed of eigenvectors of the operator s is, in fact, a common eigenstate of all 
operators Q EV . Furtheremore, the above eigenstates are orthogonal to each other: 

<k r ... k 2 k 1 \j 1 j 2 ...j r >= hi,hh 2 ,h ■ ■ ■ <5*vjr- ( 16 ) 

To find the spectrum we note that by eqs.(12) and (14) each eigenvalue of Q is 
determined only by the number of subsequent zeroes at the end of the encoding sequence 
of the corresponding eigenvector. Their multiplicities, therefore, are given by the number 
of permutations in the remaining part of the encoding sequence j\j 2 ■ ■ ■ j r - This observation 
leads to the following spectral structure of Q: 
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Proposition 1 Solution of the spectral problem for operators Q e V: 

(i) The only nondegenerate eigenvalue corresponding to the vector 1 00 . . . > has the form 

A (r) = I>; ( 17 ) 

7=0 

(ii) All eigenvectors 

\j 1 j 2 ...j r - tl ^0>, // = 0,l,...,r-l, (18) 

where j 1 ,j 2 , ■ ■ ■ ,jr-n-i belong to the set {0, 1, ... ,p - 1} and j r _ M G {1, 2, . . . ,p - 1} 
have one and the same eigenvalue 

\M = ^2 Ch . (19) 

7=0 

The multiplicity of A^ is given by mult (A (")) = (p - ly-**" 1 . 

To establish connection with the known results on the spectrum of block-hierarchical 
matrices one can use the relationship (8) between the coefficients a 7 and g 7 . By the 
application of (8) to (19) we obtain the spectrum of block-hierarchical matrices in the same 
form as in [7-9]: 



A (r) = g + (i - ^) X> 7 ^; A(0) = «o = go - gi; 



(20) 



r- 



\M = g + (l - -\ Y J P 1 q 1 -p^+i, H = 1,2, . . . ,r - 1. (21) 
3. Block-rectangular hierarchical matrices 

There are several venues to modify the block-hierarchical structure of matrices discussed 
above. The most straightforward way is based on replacing of the sub-matrices entering into 
the tensor product representation (4) of S 1 with some other matrices. For instance, one 
can substitute the identity matrices 1 with the generic diagonal matrices. This generates 
a class of block-hierarchical matrices with broken translational symmetry, where the block 
structure stays intact, while the filling of the blocks varies along the diagonal, see [12]. 
Another possibility is the replacement of matrices s. Namely, one can take a set of p x p 
matrices {sj}^ with an arbitrary content and build up a p r xp r matrix in the tensor product 
form 

(22) 
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In [13] authors introduce block-hierarchical random matrices with the blocks filled up 
randomly by zeros and ones. Each matrix from this ensemble can be generated according to 
eq.(22) if we randomly put Sj to be either identity p x p matrix or upper (lower) triangular 
matrix of l's. 

In this section we consider another modification which preserves the homogeneous block- 
hierarchical structure but turning the blocks themselves to be rectangular rather than square 
matrices as in Ex.2. 



3.1. Definition of block-rectangular hierarchical matrices 

At the first step we introduce an auxiliary translational operator T. It is instructive to 
defined it through its action on the vectors of the Hilbert space "H, see eq.(9): 

Definition 2 The translation operator T acts on vectors from H as follows 

Tuji <S> ■ ■ ■ <S> w r -i ® oj r = u r ® oj\ ® . . . ® CoV-i, 0Jj G h. (23) 



Note that in the basis (11) the operator f has the following matrix representation 



ei 



e 2 



\ 



ei 



e 2 



ei 



e 2 



V 



(24) 



ei e 2 . . . e p J 

with ei = (0, . . . 0, 1, . . . 0) T being the p-dimensional vector whose only non-zero element is 
in the i-th site while all the other elements of the matrix T are zeroes. 

Now the new family of operators V is introduced as products of block-hierarchical 
matrices and T. 

Definition 3 The family of operators V consists of all operators Q of the form 

Q = TQ, 

where Q G V and operator T is the translational operator (23). 

It is now straightforward to check that the representation of any operator from V in the 
basis (11) has a block-hierarchical structure, where blocks are rectangular matrices of the 
dimensions p 7 x p 1+1 , 7 = 1, 2, . . . , r. 
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3.2. Spectral properties of Q. 

To find the spectrum of operators from the family V we make use of a general scheme helping 
to reveal the spectrum of Q by using the properties of the operator Q . 

For the sake of convenience, first, we introduce notations for various subsets of the 
eigenvectors (15) of operator Q e V . 

(i) Let X^ denote the set of the vectors (15) having exactly \i zeros at the end of the 
encoding sequences; 

(ii) Let X^ m ' uS) be the set of the vectors (15) having m zeros in the whole encoding sequence. 
The structure of zeros is defined by the vector index v, which is the constraint partition 
of m: 

ivi = m, Vi < r — m. (25) 

i=l i=\ 

It is assumed here that zeros are gathered into t clusters, separated by non-zero elements 
of the sequence, such that the cluster of the length % enters into the sequence exactly v± 
times. Zeros standing at the end and at the beginning of the sequence are assumed to 
belong to one and the same cluster. 

The action of the operators Q and T on the above sets of eigenvectors is given by 

q X M = \(p) X M (26) 

TX^ = X^~ x \ /i^0,r; fX (r) = X (r \ (27) 

see Proposition 1 and eq.(23). From Definition 3 and eqs. (26), (27) it follows immediately 
that 

f X^x^-V, < /I < r; 
QX M = { (28) 
[ \^X {r \ ii = r. 

If ii = 0, then 

q X M = \(o) X ^'\ (29) 

where /i' is some arbitrary number less then r. 

Since the operator T possesses the property T r \ji,j 2 , . . . , j r >= \ji, J2, • • • ,j r >, each 
eigenvector of Q is simultaneously an eigenvector of Q r . In order to calculate the eigenvalues 
of <2 r , one needs, in addition, to take into account details of the encoding sequence of the 
eigenvectors: 



block-rectangular hierarchical matrices 11 

Proposition 2 All vectors x °f the set X^ m,v ^ correspond to one and the same eigenvalue 
of the operator Q r : 

Q r x = A (m; ^x, X ex^\ (30) 

The eigenvalues h.^ rn '^ are expressed in terms of in the following form: 

f (A(™))\ m = 0,r; 

A("»:") = J t (31) 

( A (o))— -Q( A W A (-D... A (i))^ ) <m<r. 

V i=l 

The corresponding multiplicities are given by 

(p - l) r , m = 0; 



mult (A< m ^) 



(f — rji — 1)1 

r(p-l) r - m 7 * ^ s n , , 0<m<r; (32) 



1, m = r. 

The particular form of the eigenvalues follows from eqs.(28), (29) with the multiplicities 
determined by the number of vectors forming the set X !<m ' v \ The latter combinatorial 
problem can be also thought as the menage problem for counting of the number of 
different ways in which it is possible to set t zeros' clusters over r — m places (m < r, 
£ = v % < r — m). The number of choices of non-zero elements in the sequence is accounted 
by the multiplier {p — i) r ^ m ; while the term r counts the number of vectors obtained by all 
possible cyclic rotations of the sequence. 

Having found the spectrum of Q we can now solve the spectral problem for the operator 
Q itself, 

= M>- (33) 

Given an eigenvector x £ X^ m '^ of Q one can construct the set of eigenvectors tj)i of the 
operator Q in the following way 



r-l 



*«= £#2 J*, 5 ; =(A("-)) 1/r e -, 1 = 1,.. .r, (34) 
\k=o J 

the £-th eigenvector corresponds to the eigenvalue A = l/5i. Note, that if the parameter r is 
not a prime number then some of the above vectors ipi might vanish. It happens when the 
primary vector x has the "period" d less then r, i.e. the following relation T d x = X takes 
place at d < r. For these vectors the formula (34) generates only d eigenvectors of Q but 
not r. To avoid a combersome formulation of the spectral theorem below we assume that 
the parameter r is a prime number. In this case the only possible values of the period d are 
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either 1 or r. Observing that two eigenvectors x an d x' generate the same set of eigenvectors 
of Q if and only if they are related by x — T k x' for some k, we obtain the following spectral 
structure of the operators from V: 

Proposition 3 Let r be a prime number, then operators Q G P have the following spectral 
data 

(i) Case m — r. The only nondegenerate eigenvalue is \( r \ It corresponds to the 
eigenvector 1 00 . . . >; 

(ii) Case m — 1,2, ... , r — 1. The eigenvalues 

A U,m;u) = ( A (m;v)^ V* . ^ ? mult ^m;u)^ = ]_ mu]t ( A (m;./)) ? j = 0, 1, . . . , 7" — 1 



correspond to the set of eigenvectors ^' m ' v ) having the fi 



orm 



r-1 _ 27ri 2(ti) 



where C is a constant fixed by the normalisation condition; 
(Hi) Case m — 0. TTie eigenvalue 

A (0,0;O) = A (0) ; ^ (A (0,0;O)) = (g ~ g ~ (g ~ ^ + p _ ^ 

corresponds to the eigenvectors obtained by application of (35) to the vectors x from the 
set X { °-'°\ 



4. Errors generation model in information sequences 

In this section we provide a schematic model of error generation in information sequences 
where the block-rectangular hierarchical matrices appear in a natural way. 

4-1. Formulation of the model 

The model is formulated in the following way. We assume that the information carrier 
("disk") has a circular shape and consists of r cells Cj, j — 1, . . . r. Each of them contains a 
digit from the set {0, . . . ,p — 1}. At the discreet moments of time t — 0, 1, . . . disc rotates 
one step clockwise. The information on the carrier can be changed between rotations of the 
disk due to an external noise which generates errors in a number of neighbouring cells. More 
specifically: 
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(a) 



(b) 




Figure 2. (a) The scheme of the model of error generation in information sequences. 
The information carrier rotating clockwise is under the influence of external noise. 
It changes the content of several consecutive cells starting from the position S. The 
level of noise defines a number of cells affected at one time-step. 

(b) The mean value of errors (solid curve) accumulated in the information 
sequence and its variance (dashed curve) after the full cycle (r = 100) plotted 
as a function of the inverse temperature f3 for the model of the Boltzmann noise. 
The data of the numeric experiment for the same model of noise are depicted by 
circles and crosses, respectively. 



- The noise affects the cells placed "rightwards" to some fixed point, S (see Fig. 2). At 
the initial moment of time t = 0, the cell Cq is located at the point S. Disk rotates such 
that at the moment of time t = £ < r the point S points to the cell Ce- 

- The intensity level of the noise determines the number 7 6 {0, . . . , r} of consecutive 
cells Ce, . . . , Ci +1 , whose information content changes randomly at the £-th time step. 
The probability to affect 7 consecutive cells is given by a 7 . The content of each affected 
cell can be replaced by any symbol from the set {0, . . . ,p — 1} with equal probability. 

In the probabilistic language each state of the disc can be described by a vector w from 
the Hilbert space % (see eq.(ll)): 



where Wj is the probability to find the disc with the content given by the sequence of 
symbols jij 2 . . . j r corresponding to j. The initial state of the system is described by the 
vector w(0) = \eji > for some j'. After a time t the system evolves into another state 





j k e{0,l,...,p-l}, (36) 



J 




(37) 
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where the coefficients Wj(t) define probabilities to find the disc having informational 
content j ± j 2 ...j r . 

The time evolution of the system can be naturally separated into two processes: rotation 
of the disc and change of the cells content due to the noise. It is easy to see that the noise 
effect on 7 consecutive cells is described by the operator 5 7 taken with the probability 
weight a 7 . The particular form of the operator s in the definition of S 1 , eq.(13), ensures 
the equidistribution of the resulting content of the affected cells. The time evolution 
corresponding to the rotation of the disc is provided by the action of the operator T. Thus, 
the result of £ steps time evolution is described by the operator 

f[f j2^A = (TQY=Q £ (38) 

i=l 7i=0 

acting on the initial vector w(0). Note, that the probabilistic interpretation of the coefficients 
a 7 's implies that Y^=o a i = 1- By eq.(17) this constraint is equivalent to = 1. 

4-2. Estimation of losses 

In this section we calculate the mean value of errors < k > accumulated in the sequence 
after one full turn of the disk and the corresponding variance < (k— < k >) 2 >. For the sake 
of simplicity of exposition the parameter p is chosen below to be equal 2. To perform the 
calculation we use the generation function Rp{a) for the moments of the number of errors 
accumulated after t steps of time evolution. The sought moments can be calculated through 
the formulae: 

< k > t=e = d a R E (a)\ a=0 ; (39) 

< (k— < k >f > t=e = d 2 a R e (a)\ a=0 - < k >\ =i . (40) 
To construct the generation function we introduce an auxiliary operator E of the form: 

E = (1 + e a P) (g) (1 + e a P) <g> . . . ® (1 + e a P); (41) 

1 = |e )(e | + |ei)(ei|, P = |e )(ei| + |ei)(e |, 

where |eo), |ei) G h and a is a parameter. In the basis (11) this operator has the following 
matrix elements 

< ej\E\e f >= e ka , (42) 

with k being the number of different symbols (compared pairwise) in the sequences j and f. 
Due to the property (42) the moment generation function can be written as 

Re(a) =<e j0 \f- e EQ e \e j0 > . (43) 
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where j encodes the initial information sequence. We apply this formula to estimate the 
mean value of errors generated after the full cycle of disk, i.e. when £ = r. Under the latter 
assumption T~ r = 1, so that 

R r {a) = 2~ r tr EQ r . 

Here we also used the fact of independence of R r (a) on the initial state jo- It is convenient 
to calculate the above trace in the basis (15) where Q and E take the diagonal form. 
Therefore, we arrive to the general expression for the moment generation function 

R r (a) = 2" r ^ A (m;,/) tr x(mi „) E = 2~ r (A (r) ) r (l + e a ) r + 2" r (A (0) ) r (l - e a )\ 



m=0 m=r 



r-1 



+2~ r J2( 1 + e a ) m (l -e a y- m J2 A{m '' u)mult ( Aim,u) ) ( 44 ) 

m=l v 

By the formulae (39), (40) taken at £ = r we obtain 

<k> = ^ (l -#(r- 1)); (45) 




r-1 \ 

< (fc- < fc >) 2 > = j ( 1 - rg(r - l) 2 + ]T «?(r - 2 - j>(j) , (46) 

j'=o / 

where we used the fact that A^ = 1 and introduced the notation g(£) = Yl i=0 

It is interesting to investigate behaviour of the moments in the case of infinitely long 
sequences. Taking the limit r — > oo in the expressions (45), (46) and assuming that 
lim^oo g(r) = exists and the sum G = YlJLo(9oo — g(j)) converges we obtain 

<k> =l(l-9oo)+O(r y t (47) 

< (k- < k >) 2 > = r - (1 - gl + 2 9oo G) + 0(r°). (48) 

Note, that both average and variance divided by the total number of cells become finite 
functions in the large-r limit. 

As an example, consider the case of Boltzmann noise where we take the probabilities 
a 7 to be proportional to e _/?7 with parameter /3 playing a role of the inverse temperature. 
From the normalising condition ^ 7 =o °7 = 1 it follows that 

1 - e-(^ _ M 

By the formula (19) the eigenvalues can be expressed explicitly as 

, , i _ e -(7+i)£ 
A,7> = L,-.^ . < 49 > 

so that 

g^ ^Wl ' 000 = 11(1-^). 

^ ' 7=1 
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This can be substituted into eqs. (47), (48). The resulting average and variance are plotted 
as functions of (3 and compared with the data of the numerical experiment on Fig.2 (b). 

One can see that at high temperatures (low j3) the mean number of errors approaches 
the half, while the variance - the quarter of the mean value of the full number of elements 
in the sequence. In this regime the noise affects almost every cell in the sequence leading to 
the Binomial distribution of the number of errors (the probability that exactly k cells change 
their content is 2 _r (£))■ This explanes the relationship between the first and the second 
moments observed at high temperatures. 

In the opposite regime of very low temperatures (large (3) the randomising events caused 
by the noise become rare and mostly affect one cell only. This results in the Poisson process 
characterised by the equality between the mean value and the variance clearly observed on 
the plot (see the right tails of the plot, Fig.2 (b)). Indeed, at low temperatures the quantity 
g(k) can be approximated by 



In the intermediate regime ((3 ~ 1) some of the randomising events start to affect 
clusters of cells rather then individual cells. This leads to some decrease of the variance in 
comparison with the one obtained from the Poisson process. An interesting feature is that at 
some point the variance reaches its maximal value and then decreases at higher temperatures. 
Qualitatively, this can be explained in the following way. At large (3 the variance grows with 
the temperature since the total number of the randomising events also increases. On the 
other hand, for high temperatures the large-scale randomising events effectively "swallow" 
the smaller ones, actually, reducing the relevant number of randomising events. Eventually, 
this process becomes dominant at very high temperatures and leads to decreasing of the 
variance. 

5. Discussion 

In the present article a new family of block-rectangular (non-symmetric) hierarchical matrices 
Q was introduced. The latter were utilised to construct the model of error generation 
in information sequences. For this model the mean number of errors and variance were 
calculated and estimated for a particular realisation of the noise. 



/3> 1 : 



goo = g(k) = l-e-P + 0(e-^), 



which, in turn, yields 
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It should be emphasised that the evolution generated by Q can be considered as a 
composition of two kinds of processes: the stochastic ultrametric diffusion described by the 
operator Q and deterministic dynamics induced by T. Since the operator T is nothing but 
the shift map on symbolic sequences, the resulting dynamics has a chaotic character [14,15]. 
In particular, it is easy to see that due to the presence of T the fully ergodig state (completely 
randomised sequence) is achieved much faster then in the pure ultrametric diffusion. The 
time necessary for the equilibration is determined by the gap in the spectrum between 1 and 
the next largest eigenvalue of the evolution operator. In the framework of our noise model 
this quantity can be estimated for the pure ultrametric diffusion as 

A (r) = 1, A (r - :) ~ 1 - (1 - e^)e~ r/3 , 

see eq. (49). The characteristic equilibration time is r cq = — 1/ ln(A^ r_1 - ) ) ~ e r/3 (l — e _/3 ) _1 
and thus becomes huge in the large-r limit. For the operator Q the closest eigenvalue to 1 
is equal to g(r — l) 1 / 7 '. When r — > oo it can be estimated as 

e ~£p < g( r - l) 1 / 7 ' < (1 - e _ ") 1/r . 

The equilibration time therefore grows linearly with r and bounded from above by f cq < re 13 . 

Generally speaking, the time evolution generated by Q. can be seen as a diffusion in 
the ultrametric space whose structure changes dynamically due to some external factors. In 
the model of error generations considered above this factor was rotation of the disk. We 
believe that similar constructions might arise in other systems having (quasi-) ultrametric 
state space, e.g. spin-glasses, macromolecules, polymers etc., subjected to time-dependent 
forces. 

Acknowledgements 

We thank Dr. S.K. Nechaev for stimulating our interest in the subject and Prof. Dr. H.-J. 
Sommers for valuable discussions. This work was supported by the Sonderforschungsbereich 
Transregio 12. 

References 

[1] G.Parisi "Infinite number of order parameters for spin-glasses", Phys. Rev. Lett. 43 1754 (1979) 

[2] M.M'ezard, G.Parisi, M.A.Virasoro "Spin-glass theory and beyond" (Singapore: World Scientific, 1987) 

[3] D. Sherrington, S.Kirkpatrick "Solvable model of a spin-glass" , Pys. Rev. Lett. 35 1792 (1975) 



block-rectangular hierarchical matrices 



18 



[4] R.Rammal, G.Toulouse, M.A.Virasoro "Ultrametricity for physicists", Rev.Mod.Phys. 58 765 (1986) 

[5] Y.Bar- Yam "Dynamics of complex system" (Wcstview Press, Boulder, Co. -1997) 

[6] G.Parisi "Complex system: a physicist's viewpoint" Physica A 263 557 (1999) 

[7] A.T.Ogielski, D.L. Stein "Dynamics on ultrametric spaces" Phys. Rev. Lett. 55 1634 (1985) 

[8] V.A.Avetisov, A.Kh.Bikulov, S.V.Kozyrev "Application of p-adic analysis to models of breaking of 

replica symmetry" J.Phys.A:Math.Gen. 32 8785 (1999) 
[9] G.Parisi, N.Sourlas "p-Adic numbers and replica symmetry breaking" Europ.Phys.J. B 14 535 (2000) 
[10] V.S.Vladimirov, I.V.Volovich, Ye.I.Zelenov "p-Adic analysis and mathematical physics" (Singapore: 

World Scientific, 1994) 

[11] A.Yu.Khrennikov, S.V.Kozyrev "Replica symmetry breaking related to a general ultrametric space I: 

Replica matrices and functionals" Physica A 359 222 (2006) 
[12] S.V.Kozyrev, V.Al.Osipov, V.A.Avetisov "Nondegenerate ultrametric diffusion" J. Math. Phys. 46 

063302 (2005) 

[13] V.A.Avetisov, A.V.Chertovich, S.K.Nechaev, O.A.Vasilyev "On scale-free and poly-scale behaviors of 

random hierarchical networks" J.Stat.Mech. P07008 (2009) 
[14] P.Pakonskil, K.Zyczkowski, M.Kus "Classical ID maps, quantum graphs and ensembles of unitary 

matrices" 3. Phys. A: Math.Gen. 34 9303 (2001) 
[15] B.Gutkin "Entropic Bounds on Scmiclassical Measures for Quantized One-Dimensional Maps" 

Comm.Math.Phys. 294 303 (2010) 



